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It was conjectured recently that Statiscally Preserved 
Structures underlie the statistical physics of turbulent trans- 
port processes. We analyze here in detail the time-dependent 
(non compact) linear operator that governs the dynamics of 
correlation functions in the case of shell models of passive 
scalar advection. The problem is generic in the sense that 
the driving velocity field is neither Gaussian nor 5-correlated 
in time. We show how to naturally discuss the dynamics in 
terms of an effective compact operator that displays "zero 
modes" which determine the anomalous scaling of the corre- 
lation functions. Since shell models have neither Lagrangian 
structure nor "shape dynamics" this example differs signifi- 
cantly from standard passive scalar advection. Nevertheless 
with the necessary modifications the generality and efficacy 
of the concept of Statistically Preserved Structures are fur- 
ther exemplified. In passing we point out a bonus of the 
present approach, in providing analytic predictions for the 
time-dependent correlation functions in decaying turbulent 
transport. 



I. INTRODUCTION 

Turbulent transport processes refer to the advection 
of a transported field 0(r, i) (scalar or vector) by a tur- 
bulent velocity field u{r,t) 0,0]. The basic equation of 
motion is linear, having the form 



9t0 = £0. 



(1) 



Here C is an operator that is built out of the turbulent 
velocity field, and as such may be stochastic. Exam- 
ples are the advection of a passive scalar 0{r,t), with the 
equation of motion 



80 



+ u-ve = KW^e 



(2) 



or a vector, like a magnetic field B{r, t) satisfying 



(9R 

— -f (m • V)B = {B ■ V)u + kV^B 



We may also consider advection, as in y], 
w whose divergence vanishes, W ■ w = 0: 



(3) 
of a vector field 



dw 

~dt 



+ {u ■ V)w = - Vp + kV'^w 



(4) 



In all these equations the velocity field u comes from 
either a solution of a fluid-mechanical equation, or is a 
random field defined with some statistical properties. A 
fundamental consequence of the linearity of the equa- 
tions of motion is that the correlation functions may be 
expressed as 



((/)(ri,i)...0(rw,t)) = 



/ni;wwpi,o)...</.(p^,o))dp 



(5) 



where (. . .) is an average over the statistics of the ini- 
tial conditions and the statistics of the advecting veloc- 
ity field. The notation r = (ri, . . . , r^) is used for sim- 
plicity. Note that we have used the passive nature of 
the transported field, i.e. the fact that the velocity is 
independent of the initial distribution of </>, to separate 
the averages over the initial conditions and the velocity. 
Such a decoupling cannot be afforded at any other time 
because of the build-up of correlations between the ad- 
vecting and advected fields. The linear operator V^, (t) 

propagates the iVth-order correlation function from time 
zero to time t. 

The evolution operator C generally includes dissipative 
terms, and without fresh input (forcing) the statistics of 
the field is time-dependent; this is the decaying case, 
Eq. (jlj). A related problem of much experimental and 
theoretical interest is forced turbulent transport where an 
input term / is added to the Eq. (^. The situations of 
interest in turbulence typically involve an input acting 
only at large scales of order L. The objects of major 
interest are the stationary correlation functions i^(^) of 
the advected field. 



FW(ri,...,rjv) = (0(ri,i)---0(rjv,t)) 



/ 



(6) 



One cares about the scaling properties at distances much 
smaller than L and at the stationary state. As usual in 
turbulent flows, the correlation functions of the advected 
field are expected to contain anomalous contributions be- 
having as 

(</)(Ari, i) • • • (/)(Arjv, t))f = A«" (0(ri, i) • • • 0(rjv, i))/, 

(7) 

with scaling exponents Cn which cannot be inferred from 
dimensional analysis. 

Recently g], two conjectures were proposed, pertain- 
ing to a wide variety of turbulent transport processes, 
without special provisos on the properties of the advect- 
ing velocity field: 

(i) In the decaying case, despite the non-stationarity of 
the statistics, there exist special functions Z^^\r) such 
that 

/W(i)= /zW(r)(^(ri,t)...(/.(r^,i))dr (8) 



are statistical integrals of motion. In the limit of in- 
finitely large system I^^^^ does not change with time. It 



follows from (^ and the conservation of /'^^ (t) that in 
the infinite size limit the Z^^^'s are left-eigenfunctions of 
the operator: 



^(A')(^) = |pW(i)^W 



{p)dp 



(9) 



note that this does not mean that the operator V ,J{t) 

admits an eigenvector decomposition, and see below for 
a further discussion of this point. 

(ii) The anomalous part of the stationary correlation 
functions in the forced problem is dominated by statis- 
tically conserved structures. In other words, at least in 
the scaling sense 



^W(r)^zW(r) . 



(10) 



A direct consequence is that the small-scale statistics of 
the transported field (p in the forced case rests on the 
understanding of the decaying problem. A by-product is 
that the scaling exponents Cn are universal, i.e. inde- 
pendent of the forcing mechanisms for any forcing that 
is statistically independent of the velocity field. 

The conjectures were exemplified in the context of shell 
models of passive scalar advection. The model's equa- 
tions read j|,^ 



d9m 



i{km+l9m+lU 



m+1 



Kkf^e^, (11) 



-'in.m' ^m' 



where the variables u„ are generated by the "Sabra" shell 
model 0] 



= i(aA;„+i-u„+2U*i+i + &fc„u„+iw*_i (12) 

+ ckn-lUn_lUn-2) - vk^Un + /„ . 



du„ 
~dt 



Here the coefficients a, 6, and c are real. In Eqs. ( p_l[ ) and 
( p^ ) the wavevectors are fc„ = /cqS". The velocity forcing 
/„ is limited to the first shell n = 0. For k — ly — 
and a + b + c ~ the energies ^„ |u„p and ^„ |6'„p 
are dynamically conserved, i.e. realization by realization. 
The statistical physics of this model were studied care- 
fully 01 in the regime of b « —0.5. Taking the forcing 
to be random (with random phases) leads to non-trivial 
statistics of the velocity field, with anomalous exponents 
that characterize the scaling behavior of the correlation 
functions. 

The operator P^^^ of Eqs. (|), (|) takes here the 
explicit form 



where n = (ni , . . . , n^) and 



.(^|0)) 



i?„,„.(i|o) = r+ 



exp[ / dsC{s) 
Jo 



(13) 



(14) 



with T"*" being the time ordering operator. Note that for 
notational simplicity we dropped the dependence on the 
initial time from V^ ' , but left it, for future purposes, in 
R{t\0). 

To demonstrate the statistical conservation laws, two 
things were done g]. First the forced problem was con- 
sidered, adding random forcing to Eq. (O): 



d9m 

IT 



C 



rn.rn' ^m' 



frr 



(fmit)f*it'))^C,nS,n.nS{t-t') . 



(15) 
(16) 



Due to phase symmetry constraints 0, there is only one 
non-zero second order correlation, but a number of dif- 
ferent higher order ones. For example, the correlation 
(6'„+2^^+i^^+i^n-i)/ is not zero. For concreteness we 
concentrated our attention to the following ones (we put 
a subscript / to stress that these are statistical averages 
in the stationary forced ensemble): 






FSU^i\Sn\'\0,nm\')f 



(17) 
(18) 
(19) 



Secondly, the decaying problem was examined, prepar- 
ing initial states 0„(i = 0) and following their evolution. 
The following objects were then computed: 



2) 



/W(i)^^(K(i)P)F(^ 



/W(i)^^(|0„(i)|>,„(i)|2)F(' 



(4) 
m 5 



(20) 
(21) 



/(«)(t) ^ Y. {K\'\9^{t)\'mt}\') F^'l, . (22) 



Figure 1 summarizes the results which are reprodcued 
from m. We show, for these three orders, (i) the time 
dependence of the nth order decaying correlation func- 
tions themselves, (ii) the time dependence of /*^^^(i). In 
panel (C) we show also for comparison the time depen- 
dence of J^^-* (t) if we replace the measured forced i^*-^-* 
by its dimensional shell dependence (i.e. the shell de- 
pendence if the Kolmogorov theory were right). We see 
that only the properly computed I^^^t) are time inde- 
pendent for times smaller than the large scale eddy turn 
over time, T^. The decay observed for times larger than 
Tl is simply due to finite size effects intervening when 
the decaying field reaches the largest scales. 

In trying to understand these results, it is very tempt- 
ing to interpret Eq. (^ as an eigenvalue equation, with 
Z^-'^' being an eigenfunction of eigenvalue 1. Unfortu- 
nately, the operator V^^' is not Hermitian, and in ad- 
dition it does not lend itself to an expansion in terms 
of eigenvectors and eigenvalues: it is not defined on a 
compact space. There are two "non-compact" directions, 
that of length scale and that of time. We thus need to 
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FIG. 1. Panel (A): time dependence of the decaying sec- 
ond order correlation functions (x), together with the time 
dependence of the statistically conserved quantities I^^' , (□). 
Equations (III) and (FL2) have been integrated with a total 
number of shells A'^ = 33. Time in the horizontal-axis is given 
in units of the large eddy turn over time Tl- Panel (B): the 
same as panel (A) but for the fourth order correlation function 
and with A'' — 25. Panel (C): the same of panel (B) but for 
the sixth order correlation function. Here we also present /'*'' 
when we replace the forced solution i^„ ^ j, with its dimen- 
sional prediction, (*) . In the simulations k = u = 5 x 10~^, 
a = 1, b = —0.4, and c — a + b. The wavevectors are 
kn = feo 2™ with n = 0, . . . , A. The smallest wavevector 
is given by ko = 0.05 while N defines the ultraviolet cut-off. 
As initial states distributions of 0n = were taken, except for 
n = 14, 15 where the field was initialized with a constant mod- 
ulus and random phases. The random forcing of the passive 
scalar was restricted to the first shell. 



learn how to take care of these before we can write down 
a proper theory. 

In the context of the passive scalar advection prob- 
lem, Eq. (g), these issues were solved elegantly in the 
framework of Lagrangian dynamics p-E2|. For the pas- 
sive scalar equation (y) the advected field is conserved 
along the trajectories of the tracer particles dr{t) = 
u{r{t),t)dt + \/2Hidj3{t), where (3{t) is a Brownian pro- 
cess. To know the scalar field at position r and time t it is 
enough to track the corresponding tracer particle back to 
its initial position p. The evolution operator V , (t) in 

Eq. (ra) coincides then with the probability density that 
N tracer particles reach the positions r at time t given 
their initial positions p. For example, to understand the 
exponent Ca one needs to focus on the dynamics of three 
tracer particles. Obviously, three particles define at any 
moment of time a triangle, which in turn is fully charac- 
terized by one length scale R (say the sum of the lengths 
of its sides), two of its internal angles, and all the an- 
gles that specify the orientation of the triangle in space. 
When the particles are advected by the turbulent veloc- 
ity field, the scale R of the triangle and its shape (angles) 
change continuously. The statement that can be made is 
that there exist distributions on the space of the triangle 
configurations, that are statistically invariant to the tur- 
bulent dynamics [p|- [l0| , |l2| . In other words, if we release 
trios of Lagrangian tracers many times into the turbulent 
fluid, and we choose the distribution of their shapes and 
sizes correctly, it will remain invariant to the turbulent 
advection [|3[ . Such statistically conserved structures are 
the aforementioned zero modes and they come to domi- 
nate the statistics of the scalar field at small scales. The 
anomalous exponents of the zero modes, such as C3, can 
be understood as the rescaling exponents characterizing 
precisely such special distributions. Of course, the same 
ideas apply to any order correlation function with the ap- 
propriate shape dynamics. The relevance of Lagrangian 
trajectories can be also demonstrated for the magnetic 
field case (g), by adding a tangent vector to the tracer 
particle, and see |jl4] for more details. 

The problem of non-compactness due to the explicit 
time dependence of the operator is taken care here by ex- 
pressing time in terms of a single scale variable R, using 
the Richardson law of turbulent diffusion [Q . Then in- 
stead of looking at the problem on the non-compact space 
of particle separation, one focuses on the space of shapes 
which is compact, and in which one can demonstrate the 
existence of eigenfucntions and eigenvalues 1 11 1^ . Ob- 
viously, for the case of the shell model considered here we 
cannot repeat verbatim the same ideology. There are no 
"shapes", and it is not immediately obvious how to re- 
late time to scales. The Lagrangian invariance is broken 
by the discretization of shell space, and the genericity of 
the time properties of the velocity field does not allow 
explicit calculations of the operator "^.^[^(i)- 

The aim of this paper is to achieve the equivalent un- 
derstanding for the shell model, that in B was orginally 



chosen to be as far removed as possible from the con- 
tinuous passive scalar problem. We will discover that 
also in this case there is a typical "moving" scale that 
carries the explicit time dependence. By considering the 
relevant operators with shell indices expressed in terms 
of the moving scale, we compactify the picture with re- 
spect to its time dependence. Moreover, in this moving 
frame we will discover that the operators decay rapidly 
as a function of shell differences. This will allow us to 
compactify the theory altogether and to offer a satisfac- 
tory understanding of the existence of the statistically 
preserved structures and its implication for the forced 
problem. 

In Sect. 2 we present the theory for 2nd order objects. 
On the basis of numerical simulations we offer an ana- 
lytic form for the operator P"-'. We show that it has an 
explicit time dependence in addition to a dependence on 
a moving scale that we identify analytically. In Sect. 3 
we use the explicit form of 7-'"' to explain why /^^^ is 
a statistical constant of the motion. The basic property 
that is crucial is the effective compactness of the opera- 
tor in the space of shells, once it is expressed in terms of 
the moving scale. Next we show how the forced station- 
ary correlation function F^'^' is obtained by solving the 
forced problem with the same propagator T^ ' . Finally 
we derive the fact that .F^^' acts as a left-eigenvector of 
P"-* with eigenvalue 1. To help in throwing light on some 
issues we also consider in this section a simple model ob- 
tained by replacing the Sabra model for the velocity field 
by a delta- function correlated field (the Kraichnan shell 
model [plJlq]). In Sects. 4 and 5 we turn to a discussion 
of the 4th order objects. We proceed in parallel to what 
had been achieved in Sects. 2 and 3. We first derive, 
on the basis of simulations and the fusion rules |1^ , the 
analytic form of V^ ' . Using this form we explain why 
/(^^ is a statistical constant of the motion when the sta- 
tionary correlation function i^'^-* is identified with Z^'^'. 
Last we turn to the forced problem, and demonstrate 
that F^"') is indeed the forced solution. This calcula- 
tion is not trivial, calling for a careful discussion of the 
time-decay and decorrelation properties of the operators 
Rn,m(t\0). Throughout the discussion we make use of the 
simpler Kraichnan shell model in which the operators are 
all computed analytically (see Appendix) to further our 
understanding of the generic case. In Sect. 6 we present 
a discussion and a summary of the paper. One very im- 
portant conclusion is that we can in fact offer an analytic 
solution for the time- dependent correlation functions in 
the decaying case; this is a considerable bonus of the 
present approach. 



II. THE FORM OF THE 2ND ORDER TIME 
PROPAGATOR 




FIG. 2. Typical time dependence of one column of the sec- 



ond order propagator 'P„|^(i)- 



Shown here is V 



(2) 



n|2oWforthe 
different times displayed in the inset in units of tq. Note that 
the maximum moves in time to lower shell numbers. 



A. Simulations 

In this section we analyze the form of the 2nd order 
propagator that governs the dynamics of the second order 
passive structure function. It is defined by: 



l^«WP>=E^SmW<l^'»(0)l' 



(23) 



The (. . .) average is over realizations of the velocity field 
and the initial conditions of the passive field. As men- 
tioned above at time t = the statistics of the advected 
field is independent of the statistics of the velocity field. 
Using simulations we can generate the matrix representa- 

(2) 

tion of 7'^|,„(i) column by column by initiating a decay- 
ing simulation (without forcing) starting with (5-function 
initial conditions in shell m. Measuring (|6'n(i)P) and av- 
eraging over many realizations of the Sabra velocity field 
we collect data for V , (t). 

In Fig. we show a typical column of V^.l^it), where 
m — 20. We used 28 shells in both Sabra and the passive 
field, with the dissipative scales being around n = 25. 
We observe two effects. First, the overall area under the 
curve decreases with time, this is the effect of the dis- 
sipation. Second, the maximum in the curve shifts to 
lower shell numbers. These are the two issues that we 
need to tackle, the time dependence and the increase in 
length scale (or, equivalently the decrease in shell num- 
ber) , which contribute to the non-compact nature of our 
operator. 

In attempting to contain these two issues we try the 
following ansatz for the propagator: 
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FIG. 3. A plot of tr^^^loit) as a function of n - m(f,20). 
The quality of the data collapse is deteriorating at the right 
tail because of viscous effects, where power law scaling crosses 
over to exponential decay. 



V. 



SiW = ^ff("-m(i,m)) 



for t > r„ 



(24) 



where r™ is a typical time scale associated with the shell 
in which the simulation was initiated. We use below 



-mC; 



/fcoTlf 



uo\ 



(25) 



with ^2 being the scaling exponent of the second order 
structure function, cf. Eq. (M). Accordingly, all times t 
below are also measured in units of tq = l/fco\/(|Mop)- 
The function H{x) has a peak at x = 0, with H{0) = 1 
and for a; > it has the form 



H{x) 



-C,2X 



(26) 



,(2) 



The location of the maximum of T-" / it) is ihit, m), and 
is a real valued function of time and of the initial peak 
location for i = 0, which is m. For t > it satisfies 
m(i, m) < m. 

To show that the ansatz (g4|) is well supported by 

the data, we show in Fig. ^ ^^„|m(*) ^^ ^ function of 
n — m(i, m). The quality of the data collapse speaks for 
itself. We draw the attention of the reader to the fact 
that the function shown in Fig. falls off sharply around 
the maximum. This will be the clue to understanding 
how to remove the noncompact dependence on the ever 
increasing scale ifi{t,m). Sums over n will be extended 
below from — c» to oo, with impunity. The important 
conclusion is that to a good approximation the following 
sum 



'H{n 



^ H{n - m) = ^ H{k) ^ Const 



(27) 



is time independent. 



B. The time dependence of the maximum 

Next we want to find an analytic expression for the 
moving scale m(t, m). In order to find the time behavior 
of the peak we examine Eq. (E3) for an initial condition 
(l^fe(O)P) — Sk^m- On the one hand, for these initial 
conditions, after time differentiating we get: 



dt 



E(i^»wi')-|Enliiw 



On the other hand, using Eq. (nAv) one finds 



^E(i^«wi^) 



2K\ki{\0^it)\' 






(28) 



(29) 



To evaluate the sum on the R.H.S of Eq. (p9), we note 
that for this linear problem the shell d from which the 
dissipation of the scalar becomes significant is indepen- 
dent of the scalar value (and thus time independent). We 
can estimate it by comparing the terms on the RHS of 
Eq. (0): 



KK^l 



Udkd 



(30) 



We can now estimate the scalar dissipation, under the 
approximation that it takes place mainly in shells with 
m> d. In this region the value of (|f?„(t)p) begins to fall 
off exponentially with fc„, and the sum in Eq. (M) is well 



approximated by the first term Kfc^(|6'rf(i)p) 
the functional form of T'"-' given by Eq. (|24 
m-m) we get: 



Plugging in 
), using Eqs. 



Y-p(2) 



dt ^-^ 



nlJt) = -J^l^Hik) 



E 

k 



CKk, 



dTrn r,-C2{d-m(t,m)) 



t 



(31) 



Examining Eq. (p^) we conclude that in order for the 
RHS to scale like t~^ for t ^ Tm, while demanding that 
for t K Tm rh{t, m) « m, we must have: 



m(i, m) = m 

m = 1 , 



l0g2 



1 

lim g(x) 



where r^ was defined in Eq. (|2£ 
we will use 



(32) 
Thus for large times 



m(t^ m) 



-;^log2( — ) 
C2 To 



t> r„ 



(33) 



Note that for large time [t ^ r^), m{t,m) becomes in- 
dependent of m. This is approporiate, since the expo- 
nential increase in typical time scales t™ when the shell 



index decreases implies that the position of the maximum 
becomes independent of its initial position. We can now 
express the time dependence of the operator "P"' in Eq. 
(P4|), solely through the time behavior of rh{t, m), by in- 
verting Eq. (p2), to find t: 



V'-V it) ex 2-^^("~™)i7(n - m) , t> 



^(2) 



(34) 



Having done so, we have gotten rid altogether of the ex- 
plicit time dependence of ■p^,^(i). Note that the depn- 
dence of the operator on both its shell indices turns nat- 
urally to a dependence on the difference between these 
indices and the single moving shellThis is the first impor- 
tant step in overcoming the non-compactness of our 2nd 
order operator. 



III. CONSEQUENCES OF THE FORM OF THE 
2ND ORDER PROPAGATOR 

At this point we can reap the benefit of the explicit 
form of the 2nd order propagator Eq. (p4). First we de- 
rive the existence of the statistical constant of the motion 

/(2). 



A. Second order constant of the motion 

Returning to the definition of /'^^\ Eq. (EQ), and rec- 
ognizing that Fn (X 2~^2" (which is also demonstrated 
in the next Subsection), we see that we need to evaluate 
the weighted sum E„(l^n(0P)2"^'"- Since the problem 
is linear, any initial condition can be represented as a 
weighted sum of (5- function initial conditions, and there- 
fore we only need to consider sums of the form 



■> 

it 
I 

c 

fN 
I 

'^ 



y p^i^ it)2- 



<2n 



2-C2m Y^ jji^^ _ jj^^2' 



C2(n- 



(35) 



As the components of the sum are a function of n only 
through the combination n — ih, we can change the sum- 
mation to run on n — rh. In light of Eq. ( p7| ) the sum is 
time independent. In Fig. |j we show the summand as a 
function of time and n — fn. 



B. The Forced 2nd order steady state solution 

For the forced solution we can use again the fact that 
the statistics of the velocity field has no correlation with 
the forcing of the passive scalar field at any time. There- 
fore we have: 



{\On{t)\' 



•^0 J, 



t')(|/fc(i')P>rfi' • (36) 
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FIG. 4. The summand H{n - m)2~'^2("-™' as a function 
oi n — m. 



We should think of this equation only in the limit of 
t —> oo, since we need to eliminate the effects of expo- 
nentially decaying initial value terms that do not con- 
tribute to the stationary forced correlation function. 
With a force which is Gaussian white noise, we write 
(l/fe(i)H = Ph,m- Using Eq. (H) for the propagator 
we get: 



(l^nWP)/«/' 



H{n-m{t-t',m)) (37) 



X 2-C2(m-m(t-t\m)) ^^1 



We remind the reader that r^ is the time it takes for the 
initial (5-function to develop a "scaling tail" for n > m, 
and now m is the shell at which the random forcing is 
localized. The idea here is to use the fact that we know 
how to eliminate the time variable in favor of the moving 
scale variable rh(t, m). Changing variables of integration 
to m, using Eq. (p2) we can write explicitly for rh <^ n: 

(|f?„(t)n/cxC2ln(2)/2 (38) 

2-C2{n-m) 2-C2{m-m) 2-C2rh ^^^ 



Note that we have extended in a formal manner the range 
of shell indices all the way to — cxd, to allow for a long 
development of a self similar solution. Naturally, since 
the integral converges quickly, this is immaterial. Finally, 
using (O): 



Fi^) ^ Const X 2-^^^" . 



(39) 



This solution has the expected 2 ''^", and is time inde- 
pendent. 



We note at this point that the forced solution Fn ' 
had been shown to be a left eigenfunction of eigenvalue 
1 in Eq. (|3q). Thus the first two subsections together 
fully demonstrate the two conjectures (i) and (ii) from 
the introduction for the case of the second order objects. 



C. Why this simple time dependence? 

The knowledgeable reader might have noticed at this 
point that the explicit time dependence of the 2nd order 
propagator, as displayed in Eq. (gj) is very simple. The 
exponent of time, i~^, is not anomalous, and appears 
independent of the second order exponent of the velocity 
field. This is not so in the understood example of the 
Kraichnan model of passive scalar advection, in which 
the time dependence of the operator is anomalous ||l5|,|n| . 
To clarify this point we turn to the analysis of the passive 
scalar shell model driven by a (5-correlated velocity field 
[g[ . In other words, for the velocity field u in Eq. (pT|), we 
use a Gaussian field, (5-correlated in time, that satisfies: 



(u„(tX(t')>='5„,r„^(i-t')C„ 



(40) 



The calculations are described in Appendix |^, with the 
following results: 



-(i0„(i)n=A42)^(i0„(i)p) 

Where the matrix M^'^^ is given by 



Mi 
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^n ^ T-„+l 



Tn+1 



(41) 



(42) 



Here t„ = 2 (^ ^^"//coCo. Since this matrix is time in- 
dependent we have 



,(2) 



r::^:^{t) = [eMtM^n 



(43) 



It is straigthforward to check that the 2nd order forced 
solution (1 6'„ 1 2)^ - 2-(2-«)» ~ r„ is a zero mode of M^^'. 
In this case it is also straighforward to prove that /^^^ in 
Eq. ( poj ) is a conserved variable (in the infinite system 
limit). To do so we note that on the one hand from Eq. 
([A5|) we have the following exact equation: 



dt ^ Trf+i 

n—l 



(44) 



The explicit form of the quantity J*^^^ is in this case 

/^'^=E^™(I^™WI') (45) 

m 

The rate of change of this object is 



;^E^"(I^™WI') -T.^mM^^UK{t)\') (46) 



dt 

m m,n 

Using the properties of M^ ' we can write: 
d 



dt 



Y.r,M^{t)?) = -{\Oa{t)? 



(47) 



Taking the ratio of Eq. (||) and Eq. (|^ we see that 
for the limit d — *■ cx) the quantity I^"^' is conserved with 
respect to the sum X]m(l^™(^)P)- 

Now write the propagator in the form pertianing to 
t > r™. 



V^ 



(2) 



Xi)=c(^)"i/(n-m(t,m*)) 



with 



m(i) = mo - ^--- log2(t/T,„) 



(48) 



(49) 



Write now the conservation law just proven as 

I^^^=T.^nVl^in^t)^Gor.Si. (50) 

n 

Using the form (^8|) we require 

E H{n - m(t, TO*))2-(2-«)(«-"™(*.™*)) = Const . (51) 

n 

Obviously, this is constant iff a = 1, demonstrating the 
point that the explicit time dependence in our propagator 
is not anomalous. 

IV. THE 4TH ORDER PROPAGATOR 
A. Simulations 

The 4th order propagator is defined by: 

(ien(i)n^™(i)p) =pS.i,,,(o(i^p(o)ne,(o)p) . (52) 

We remind the reader that the LHS has also contribu- 
tions from other initial con- 
ditions, i.e. (6'„+2(0)6i;+i(0)6i;+i(0)6i„_i(0)) but these 
contributions appear in the numerics to be very small, 
and will not be considered in this paper. For (S-function 
initial conditions (say on shell p) it is sufficient to con- 
sider V I (i). For ra,n ^ n.q and for large times, 
V i it) IS indistinguishable from V , it). 

n,m\p,q^ i '^ n,m\p,p^ ' 

First we studied the typical time dependence of the 
operator via direct simulations. In Fig. g we plot the 
diagonal elements V\ ^35 25 (^) ^^ functions of n for dif- 
ferent times. Note the movement of the peak and the 
decay of the function. This is very similar to what we 



,xlO 




,xlO 




FIG. 5. The diagonal elements of P^ ^j^ 25(*) ^^ functions 
of n for five different times. The simulations were performed 
with 30 shells. 



20 -10 __0 10 20 

n-m(t,25) 

FIG. 6. The diagonal elements of t'^''^'^^Pi'^i|25,25 
function of n — rh{t, 25). 



it) 



found for the second order propagator. In order to pro- 
ceed we need to guess an analytic form for the propagator 
and compare it with the numerical data. 

Our ansatz for the 4th order propagator is constructed 
using the fusion rules llfl]. For the forced 4th order cor- 
relation functions the fusion rules predict that asymptot- 
icaly for |n — r7?,| ^ 1 



^(4) ^ 2"C4min(m,n)2~C2|m-n 



(53) 



This form was amply tested and demonstrated for shell 
models in jlj] . It was shown that the asymptotic form is 
obtained very rapidly, for any |n — to| > 1. Accordingly 
we expect that 






C4/C2 



G\ min(?n, n) — m(t,p))2 



"C2 |m— n 



(54) 



where the function rh{t,p) is the same as in Eq. (pi) but 
with p replacing m. The function G{x) is expected to 
have, for a: ^ 0, the scaling form: 



G{x) ex 2"'^''^ 



(55) 



The form (p3) is very well supported by the data. In 
Fig. y we replot the data of Fig. || multiplied by t'^'>-/i^ , 
as a function of n — m(i, 25), where m(i, 25) solves Eq. 
(p^). It is obvious that the form ( |54|) is justified for the 
diagonal. 

It is more difficult to demonstrate the full tensor by 
direct simulations; the off-diagonal elements are more 
noisy, and the scaling behavior is somewhat less apparent 
than on the diagonal. We can however obtain much bet- 



ter data for the Kraichnan model, for which V 



(4) 



T-\P,P 



it) 



can be computed essentially analytically. In Appendix 
H we present the derivation. Here we show in Fig. M 
[t] for three different times. The spread and 



V. 



(4) 
n,m|18,18 



decay are apparent. In Fig. |^ the same data is shown 
after multiplying it by i''*/'»^, and replotting it as a func- 
tion of (n — rh{t, 18), m — rh{t, 18)). Now the function is 
preserved with respect to time. 



V. CONSEQUENCES OF THE FORM OF THE 
4TH ORDER PROPAGATOR 

A. The 4th order constant of the motion 

According to the conjectures discussed in the intro- 
duction (in particular Eq. (|lO|)), we expect the forced 
solution _F(^^ to act as the left eigenfunction of eigen- 
value 1, Z'^. Here we demonstrate that /^^^ as defined 
by Eq. ( pl| ) is indeed a constant of the motion. Using 
for F^^^f Eq. (53), and expressing t^^/'^^ jjj terms of rh we 
get: 



/'^ ' (t) — 2^ GI min(TO, n) — rh{p, t) 

^ 2^2C2|ni-n|2-C4(min(m,n)-m(p,t)) (5Q) 



As in Eq. ( p5D , the time dependance of the sum is elim- 
inated because time is introduced only through the ex- 
pression min(?n, n) — rh{p,t). Consequently the object 
/(^^ becomes time independent. We demonstrate this in- 
variancc for the diagonal part of the summand in Fig. 
H. To display the invariance for the whole weighted ten- 
sor we employ again the data presented in Fig. 0. After 




,(4) 



FIG. 7. The logarithm of the elements of P^^^Mg isi^) ^^^ 
the Kraichnan shell model as a function of n and m for the 
three different times 1.54 x 10~''ro (panel a), 1.67 x 10~^to 
(panel b) and 1.74 x lO^^ro (panel c). 




FIG. 8. The logarithm of the elements 



of 



^Ci/<2-py^ )^ 18 (i) for the Kraichnan shell model as a func- 
tion of (n — m{t, 18), m — ih{t, 18)) . The times are the same 
as in Fig. M. The invariance of the function is obvious. 
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FIG. 9. The weighted elements Fn^nV^ ^|-^g 15 (^) ^^ ^ func- 



tion of 71 — m{t, 15). 



Ti,n| 15,: 



(4) 



multiplication by the weights Fn,m and replotting in mov- 
ing coordinates, the constancy of the summand of /'^' is 
demonstrated. This is done in Fig. nO, using the analytic 
results of Appendix ^. 



B. The forced 4th order steady state solution 



Finally, we can calculate the analog of Eq. (pq), for 
the steady state 4-point function in a system forced by 
gaussian white noise. Returning to Eq. (14) we write 



i^(4) 



dsi ■■ ■ ds4{Rn,pit\si)R^p,{t\s2) 



i?„,,(t|s3)i?;,,,,(t|54))(/p(5i)/;(52)/,(s3)/;,(s4)) , (57) 

where we have used the statistical independence of the 
forcing from the velocity field. We note that in Eq. ( |57| ) 
the time integration can be (and should be) extended to 
arbitrarily long times, to get a stationary forced correla- 
tion function. This way we also get rid of exponentially 
decaying initial value terms. Using the correlation prop- 
erties of the forcing Eq. ( |lq ) we obtain 



Fi%it)=2CpCg [dsi [dS2{\RnAt\si)\''\RmAt\s2)\' 
JQ Jo 



(58) 

We next split the integral into two domains in which 
S2 < si and the opposite. Consider the first domain in 
which the integral on the RHS has the form 

[dsi [ dS2{\RnAt\si)f\Rm.At\s2)f) 
JQ JO 

= Ll [dS2{\Rn,p{t\si)\''\RmAt\si)\''\RiAsi\s2)\'') ■ (59) 
JO JO 




FIG. 10. The elements of -F'i*m7'i'^L|i8,i8(*) ^°^ ^^'^ Kraich- 
nan shell model as a function of (n — m{t, 18), m — rh{t, 18)) . 
The times are the same as in Fig. m The invariance of the 
function is obvious. 
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To proceed we need to consider the decay time of the 
operator i?„,m(^|io) compared to the decorrelation prop- 
erties of products of such operators at different times. 
On the one hand wc know that these operators depend 
exphcitly on time, decaying hke a power of time (cf. Eq. 
(p4|)). On the other hand, we expect the correlation of 
products of different time operators to decay exponen- 
tiaUy, since the operators Rn.,m{t\to) contain the chaotic 
velocity field that appears in the exponential, cf. Eq. 
(|lj). The time domain is arbitrarily long, but through- 
out most of the time integration the product is actually 
decorrelated, and we can write 



dSi / dS2{\Rn,p{t\si)?\R^AAsi)?\Ri,q{si\s2f 



dsi{\Rn,p[Asi)\^\RniAAsi)?) / ds2{\Ri,q{si\s2)\ 



We can now perform the integral over S2, with a result 
independent of either si or t. Finally, if we choose the 
forcing in Eq. (5^) as a single shell forcing on the shell 
P, {\fk{t)?\fi{ty) ^ C^5k,p6i,p we get: 






-C2I™— n 






C4/C2 



(61) 



where we have used the analytic asymptotic form of 



V, 



(4) 



I (i), Eq. (|54|). Changing the integration variable 
to m{t — t',p) we get: 






C2 l^^~^l c\~ C,i^{yii\v^{ra ^ny) 

22C4A2-C2m^^ , (62) 

00 

or, we find a time independent solution: 

F^d^ = Const X 2-C2|™-«|2-C4mm(m.„) (g3) 

As expected, this is the correct form of the 4th order 
correlation function, in agreement with the fusion rules. 
The theory for the sixth and higher order correlation 
functions follows the same lines, ane will not be repro- 
duced here. 



VI. SUMMARY AND CONCLUDING REMARKS 

In summary, we examined in detail the statistical 
physics of the shell model of a passive scalar advected 
by a turbulent velocity field. We presented a theory to 
explain and solidify the two conjectures proposed in |^] 
and reproduced in the introduction. These conjectures 
state that (i) in the decaying problem there exist in- 
finitely many statistically conserved quantities, denoted 
above as I^^^^; (ii) These quantities are obtained by in- 
tegrating (or summing) the decaying correlation func- 
tions against the stationary correlation functions of the 



forced problem. We have pointed out that the conjec- 
tures imply that the forced solutions are /e/t-eigenvectors 
of eigenvalue 1 of the propagators "P^ '. For the model 
discussed above we have established these conjectures by 
examining the form of the propagators P*- ■*. Using nu- 
merical simulations as a clue, we proposed analytic ex- 
pressions for the operators P'^' and V^ \ pointing out 
that similar concepts (fusion rules |16[| in particular) can 
be used to write down also the higher order operators. 
We checked the analytic forms against the simulations, 
and proceeded to demonstrate that the forced, station- 
ary correlation functions arc indeed left-eigenvectors of 
eigenvalue 1 of these operators. This implies that the 
objects 7*^^^ are indeed constants of the motion. Next 
we derived the forced, stationary correlation functions, 
/„„and showed that the form of our operators dictates scal- 
ing solutions in agreement with the fusion rules. As a 
result the two conjectures were confirmed. In our ana- 
lytic calculations we used repeatedly the fact that the 
operators "compactify" in shell-space once expressed in 
terms of a single moving scale whose dynamics was de- 
termined analytically. 

One should state a caveat at this point: the analytic 
form of the opeartors 'P and 'P ' was guessed on the 
basis of numerics and the fusion rules. Although they 
appear to agree with the simulations, we cannot state 
that the forms are exact. Accordingly, until these forms 
are derived from first principles, the exact status of the 
conjectures is not established. It may be that the con- 
jectures are only satisfied to a good approximation. This 
question needs to be addressed in future research. 

Notwithstanding this caveat, we should point out a 
surprising bonus of the approach discussed in this pa- 
per: we have at hand an analytic form of the propaga- 
tors. We can thus provide analytic predictions for the 
decaying correlation functions for arbitrary initial condi- 
tions. Considering that the velocity field is a solution of 
a highly non-trivial chaotic dynamical system, and that 
the passive scalar is slaved to it, it is quite gratifying 
that nevertheless one can offer analytic solutions for the 
time-dependent correlation functions of the latter. It is 
of course very tempting to hope that a similar theory 
can be developed in other cases of turbulent transport, 
leading to analytic predictability of the time-dependent 
correlation functions in the decaying case. Since this pa- 
per demonstrated that the Lagrangian structure is not a 
prerequisite for the exisitence of Statistically Preserved 
Structures, we feel that such a theory should be sought in 
the Eulerian frame in which calculations are much easier 
than in the Lagrangian frame. This development should 
be addressed by future research. 
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— (|e„(t)P) = C„+ifcf,+i(|0„+i(t)|^> +C„A:J;(l^n-i(t)|'> 
-(C„fc2+C„+ifc^+i)(|0„Wn . (A5) 

This can be written in matrix form as: 



-(|&„(i)P)=Af(2)^(|0„(t)P) 



(A6) 



Where the matrix M^"^' is given by Eq. (M2). It is time 
independent and thus a solution for ^„|„(i), defined in 
Eq. @, can be written as Eq. @ 



APPENDIX A: THE KRAICHNAN SHELL 
MODEL 



2. The 4th order operator 



For the velocity field u in Eq. ( ]11| ) we use a Gaussian, 
delta-correlated in time field that satisfies: 



{Un{t)u*„^{t')) = 5n.,yn5[t - t')Cn , 



(Al) 



For this simple model we can find a closed form equation 
for the time derivative of T'^'^\t). For simplicity we set 
the diffusivity k = 0, and replace its effect by truncating 
the opeartor at the dissipative shell d (cf. Eq. (|3C|)). 



1. The 2nd order operator 

We evaluate the 2nd order propagator's time derivative 
by multiplying Eq. ( |ll[ ) by 0* and adding the complex 
conjugate to get: 



-(|0„(i)P) = *A;„+iK+i(i)0„+i(t)C(i)) 

+iK{ui(t)ei{t)en-i{t)) + C.C , 



(A2) 



Using Gaussian integration by parts we compute the 3rd 
order correlation functions including the velocity: 

{ei{t)Om{t)Ura{t)) = ( dt' Y,{u^{t)u;{t')) 



(A3) 






From Eq. (O) we have for the functional derivatives: 
SOpit) 



5ul{t') 
Su*{t') 



ie{t~t')6p,gkgeg_l{t'), 
-iQ{t~t')6p+l,gkg9l{t'), 



(A4) 



where 8(t) is the step function, 8(t) — 0,t < 0, 0(i) = 1, 
t > 0, 6(01= 1/2. Plugging Eqs.Q, (|A§ into Eq. 
(^), Eq. (|A2|) becomes 



Let us consider the propagator of the 4-point correla- 
tion function (|6'„(i)|2|6l„(t)|2) : 



dt 



{\0,M'\O.nitr) = Mi%^^^^{\9,it)\^\9g{t)\^) , (A7) 



where the operator M^ ' can be computed in analogy to 
Eqs. (PI), (ph, (11): 



n,rn,'p,q 



(2) 



' n 1 ^n.m '^n.m,-{-l 1 1 '^n.p'^n — l.q i ^n,q'^n — l.p 1 
+ '^'''n+l\^n,m—5n+l.mj{5n,p5n+l,q+5n,qSn+l,pj ■ (A8) 

We note that M^'^' is not symmetric under the exchange 



of left and right indices, i. e. Mn,in,p,q ^ Mp, 



and 



thus admits different left and right eigenvectors. The 
zero-mode of Eq. ( A7) satisfies 



m(4) y(4) = 

n,7n.p,q p,q ^ 



(A9) 



and is expected to be a symmetric function of the form 

y1% = 2-'^^"""('"'")/^d"^ - "I) ■ (AlO) 

Equivalently one can consider a left zero-mode of M^ ' , 
which we denoted above by Z^^^ (cf. Eq. (||)), 



m(^) Z^'^^ = 

p,q,n,m p,q 



Z(4)^ = 2-«- --("■") /L(|m-n|) 



(All) 
(A12) 



We will show that both left and right zero modes have 
overall scaling exponent ^4, multiplied by a function 
/^/L(|to - n|), which scales like 2-?^l™-"l provided 
Itti — n| 3> 1, in agreement with the fusion rules. We 
therefore propose the following ansatz 






(A13) 
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Plugging this ansatz into Eqs. ([A9|), (All), we find three 
different cases : (i) to = n, (ii) to, = n ± 1 and (iii) 
|m — n| > 1. This last case, which is identical for both 
left and right equations, reads (assuming tti > ?i + 1) 



[6] 

[7] 

[8] 
[9] 

[10] 
+T-l^2-^"f^/'^{m -n-l) + T-^2'^^f^/'^{m - n + 1) , (AM) 

which, defining P — (4 — 2^2, yields the following recur- 
sion relation for the coefficients a ■ in Eq. (A13), >- ' 



irn'+' 



'ri+1 



'm+ 



^)f^/^{m-n) 



'rn+ 



i/«/L(TO-n+l) + r„V^/^(r 



R/L 



2-^' 



3-2 



^1 ' - rj-i - r, 1 



2^'0-3 R/L . . „ 



"j-i 



(A15) 



It then remains to determine /3, /^^'"(O) and Qi , which 
is done with the help of cases (i) and (ii) above. In the 
case of the right zero-mode, we have 

(1 + T-')f^{0) - 2t~\1 + 20T-')f^il) , (A16) 
(1 + At-' + r-i)/i^(l) ^ {r-' + 2-^Ti)/«(0) 

+ (1 + 2'^)r- V"(2) , (A17) 

whereas the left zero-mode yields 

(1 + rri)/^0) = rri(l + 2''rri)/L(l) , (A18) 
(1 + ir,-^ + r2-i)/L(l) ^ 2(rri + 2-^T,)f^{0) 

+ {l + 2^)T^^f^{2). (A19) 

We note that, provided we impose f^{0) = 1/2/^(0), 
which amounts to fixing the arbitrary relative multiplica- 
tive factor between Y^ ^ and Z^ \ we obtain two iden- 
tical solutions, i. e. a^ = a^ Vj > 1. The anomaly /3 is 
then the same for both systems of equations and can be 
determined numerically. 
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